Ignition of superconducting vortices by acoustic standing waves 
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Nucleation of vortices in a superconductor below the first critical field can be assisted by transverse 
sound in the GHz frequency range. Vortices will enter and exist the superconductor at the frequency 
of the sound. We compute the threshold parameters of the sound and show that this effect is within 
experimental reach. 
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A superconducting cylinder rotated at an angular ve- 
locity about its symmetry axis develops a magnetic 
moment M = — (mc/27re)ll, where m and e are bare elec- 
tron mass and charge, and c is the speed of light. This ef- 
fect predicted by London [l| has been subsequently tested 
in experiment and proved with an accuracy to many sig- 
nificant figures. It is a consequence of a more general gy- 
romagnetic effect predicted by Barnett [1] : "A substance 
which is magnetic must become magnetized by a sort of 
molecular gyroscopic motion on receiving an angular ve- 
locity" . Barnet effect is, in its turn, a consequence of the 
Larmor theorem 0]: In the rotating frame of reference 
the action of the rotation on charged particles is equiv- 
alent to the action of the magnetic field Ho = 0/7 Q, 
where 7 is the gyromagnetic ratio. For electron's orbital 
motion 7 = e/(2mc) « 0.9 x lO"^ (Gauss) s"^ Thus, 
in practice, the fictitious field in the reference frame of a 
rotating macroscopic cylinder can hardly exceed a frac- 
tion of a milligauss. This would be well below the lower 
critical field Hd when the temperature of the supercon- 
ductor is not too close to Tc. Due to the Meissner ef- 
fect Q (considered in the frame of the rotating cylin- 
der) such a field would be expelled from the bulk of the 
cylinder by a superconducting current induced at the sur- 
face. Writing B = Ho + 47rM = for the total field in 
the bulk, one obtains the London's magnetic moment, 
M = — Hn/(47r) = —{mc/2TTe)fl. Due to the symmetry 
of the problem it is the same in the rotating and labora- 
tory frames. 

In this Letter we would like to take this problem a little 
further and look at the consequence of an angular veloc- 
ity well beyond the experimental limit. In particular, we 
are interested in the rotational velocity of a magnitude 
that would generate a fictitious magnetic field that ex- 
ceeds Hci. If Larmor 's theorem still holds, than it must 
be the case that a superconducting vortex enters the bulk 
of the cylinder. This would require the angular velocity 
to be of order 10^s~^, clearly surpassing the feasible ex- 
perimental value for a mechanical rotation. While this 
scenario is merely a thought experiment we will use it as 
a motivation to study the effect of local rotations gener- 
ated in a superconductor by high frequency ultrasound. 



Interaction of sound with vortices has been studied in 
the past S, JB, Q ■ Radiation of phonons by supersonic 
vortices Isl |9|, phonon contribution to the vortex mass 
IC , I 111 . Il2l | , and decoherence of flux qubits by phonons 
ISlllJ] have been investigated. In this Letter we are ad- 
dressing a completely different problem - possibility of 
the nucleation of a vortex by sound. 

Within continuous elastic theory, local deformations 
are described by the displacement vector field, u(r,t). 
We will be interested in the effect of transverse sound 
waves. Such waves create shear deformations of the crys- 
tal lattice, such that 



0. 



(1) 



In the long-wave limit they do not affect the density of the 
ionic lattice but result in a local rotation at an angular 



velocity [15 1 



rj(r,t) = -Vxu. 



(2) 



The frequency of ultrasound achievable in experiment 
with, e.g., surface acoustic waves can easily be in the ball- 
park of / ~ 10^°s~^ [l6|. According to Eq. ^ a sound 
of such frequency and amplitude of a few nanometers can 
provide 17 ~ 10^s~^ that can generate fictitious magnetic 
fields above Hd. For practical purposes, it may be con- 
venient to loosen the restriction on the frequency and 
amplitude of ultrasound by applying an external mag- 
netic field H near, but less than, Hd. We shall see that 
within one percent of Hd , vortices can be ignited by the 
ultrasound in the GHz frequency range. 

For a vortex to enter a superconductor, the Gibbs free 
energy of the system must be lowered. We compute the 
extra free energy due to the vortex and determine the 
condition at which it becomes negative. It should be 
noted that the system under consideration is dynami- 
cal, and therefore is not at a thermodynamical equilib- 
rium. However, we are interested in the free energy of the 
Cooper pairs which can adjust to the changes of state in 
a time scale orders of magnitude shorter than the pe- 
riod of the sound. This time scale is proportional to the 



2 



relaxation time t of the cooper pairs, i.e. t ^ 10~^^s. 
As mentioned before, the period of the sound T — 1/f 
will be always greater than 10~^°s. Under these condi- 
tions, our system is adiabatic and the thermodynamic 
equilibrium can be safely established. The calculation 
that follows is similar to the conventional calculation of 
Hci. The presence of il, however, introduces a new fea- 
ture into this calculation so we will follow it all the way 
through to show how the sound enters the problem. 

It is convenient to calculate the extra free energy in 
terms of the magnetic field and its spacial derivatives. 
The electric field produced by the time derivatives will be 
neglected. The kinetic energy of the superfluid is irisTOV^ 
where Us is the number density of the superconducting 
electrons and 



V = I — Vi^ - A 

m*c \ e* 



(3) 



is the velocity of the cooper pairs with ip and A being the 
phase of the superfluid wavefunction and the magnetic 
vector potential, respectively. The stared quantities rep- 
resent the effective mass and charge of cooper pairs. We 
will take them to be m* — 2m and e* — 2e. The normal 
electrons experience viscous forces as they move relative 
to the nuclei contributing zero average normal current. 
The ionic charge per unit volume consisting of the nuclei 
and the normal electrons is therefore exactly opposite to 
that of the cooper pairs. The total current is then 



in which po is the combined mass density of ions and 
normal electrons, Xikim is the tensor of elastic coefficients 
and Uik = ^{diUk + dkUi) is the strain tensor. Using 
Maxwell's equation V x B = {Air/cjj and combining Eqs. 
([7]) and ^ the expression for the total Gibbs free energy 
yields 



G = 



1 

8^ 
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(9) 



j = ens(v - u), 



(4) 



Here, J-q is the free energy in the absence of currents, 
fields, and sound, A = \J mc?- j Aurise^ is the London pen- 
etration depth, /(r) = ([V'l/IV'oo |)^ in which i\} is the 
complex order parameter and |'0oo| — \J ng/l is the or- 
der parameter in the absence of gradients and fields, and 
P — Po+rism is the total mass density of the superconduc- 
tor. The fourth term can be recognized as the interaction 
of the external magnetic field with the magnetization. It 
is this term that is responsible for the nucleation of vor- 
tices in the absence of sound when H > Hd- 

Before we can calculate the free energy of Eq. ([9]) 
we must first work out the magnetic field. This can be 
done by replacing the current in the Maxwell's equation 
V X B = (47r/c)j with Eq. (O and defining a gauge 
invariant vector potential Q = A — (?ic/2e) so that 
we obtain the following equation: 



where u is the velocity of ions. Eq. ^ reflects the fact 
that the electric current corresponds to the motion of 
electrons relative to ions. It is invariant with respect to 
the motion of the reference frame. With Eq. ^ in mind 
we can write the gauge invariant current in terms of ip 
and u as 



where 



riseh ( 2e . 



Ae// = Ah u 

e 



(5) 



(6) 



is the effective vector potential felt by the electrons in 
the rotating frame of the ions [13] . In terms of the to- 
tal current j, the kinetic energy of the superconducting 
electrons may be expressed in the form 



(7) 
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The energy of the sound is 
1 



Es ^ I d^r^^poii-' - XtkimUikUim) (8) 



A^V X (V X Q) + /(r)Q 



mc 



/(r)u. (10) 



For V(p = (Q = A) Eq. PH)) becomes equivalent to the 
London's equation with a source. When a vortex enters 
a superconductor the phase must be quantized accord- 
ing to the condition f V(p ■ dl = 2tt. For certainty we 
consider a transverse standing sound wave having one 
node at the center of a superconducting slab of thick- 
ness d large compared to the coherence length ^. The 
external field is applied parallel to the slab, see Fig. 1. 
In this case A^ = 2d. Generalization to standing waves 
with many nodes is straightforward. If the field is close 
to Hci, a vortex will periodically enter and exit the slab. 
The boundary condition on the current is J_l ■ n = 0, 
where n is the direction of the surface. If the thickness 
of the slab d is of order or less then A, this boundary will 
distort the cylindrical symmetry of the vortex. We can 
satisfy the boundary condition by placing image vortices 
of alternating sign a distance d apart on the outside of the 
slab. The equation for the magnetic field, in the region 
r > ^ where = 1, can then be written in two parts, 
namely B = Bq + B^, , such that the first term satisfies 



A^V X (V X Bo) + Bo 



2mc , 



(11) 
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FIG. 1: Standing wave in a slab with one node at the center. 
The vortex is generated at the node where is maximum. 



while the second is a solution of 

oo 

A^V X (V X B,) + B, = $oe. ^ (-l)"5(r + nde,), 



n— — OO 



(12) 

where $o = hc/2e is the flux quantum. Notice that Eqs. 
PT|) and (fT^ can be obtained by taking a curl of Eq. 

with the account of the vortex cores represented by 
the delta functions. 

Since we are interested in standing sound waves we can 
choose the displacement vector u to be 



u(r, t) = uo su\{kx) sin{ijjt)ey. 



(13) 



The quantity k = lo / v = 27r/As — Tr/d is the wave num- 
ber with As and v being the wavelength and the speed of 
sound respectively. It is easy to see from Eq. that Vl is 
maximum at the nodes. The corresponding solutions to 
Eqs. pT|) and with the boundary condition B = H 
at a; = ±c?/2 are 



Bo(x) 
B.(r) 
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B, 



(-l)"b(r„) 



where 
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B, 



b(r„) 



sinh(d/2A) 
sinh(d/ A) 
2mc Jl 



cosh 



e 1 
$0 



Xo(|r + ndej,\/X)ez 



(14) 
(15) 

(16) 
(17) 
(18) 



27rA2' 

and K{) is a zeroth-order Hankel function of imaginary 
argument. The first term in Eq. (|14p is the Meissner 
field while the second is due to the sound. 

Let us now integrate by parts the third term in Eq. ([9|) 
and insert B = Bq + B,j. By doing so we obtain 



G — Go 



47r 



dr'^B,; 



2mc 



rj + Bs 
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A2 



f{r. 



■(V X B„) • (V X Bs 



(19) 



where Go is the Gibbs free energy without a vortex and 

A2 



AE = ^ I (fj 



1 

An 
1 



fir) 



(V X Bm) • (V X B,) 



B 



" fir) 



(V X B„)2 



/ dr'U ■ [B, + BAfj 



(20) 



is the vortex energy. One can simplify the volume inte- 
grals in Eq. p9| by separating the integration over the 
core from the integration over the volume outside of the 
core. When the latter is integrated by parts, the integrals 
outside the core cancel, and the free energy in Eq. (fT9|l 
with the help of Eq. (fTT|) becomes 



where 



AJ^ = ATi + AJS + AT3 + AE , 



ATi = — / dr'^B,, ■ n + B, 

At: [ e 

a2 r 

AT2 ^ — * B.„ X (V X Bs) ■ ds 
47r 



AT. = 



1 



d^r 



A2 
fir) 



■(V X B„) • (V X Bs). 



(21) 

(22) 
(23) 

(24) 



The subscript c indicates an integration over the core. 
The surface integral in Eq. ([^5]) is over the boundary of 
the normal core. Near the vortex core f{r) = (r/a)'^, 
where a w ^. It is straightforward to check that in the 
limit r — > the exact solution to Eq. (fTUl) for the vector 
potential As(r) is 



A,(r) 



2mc 



(25) 



Then the magnetic field Bg = V x As(r) generated by 
the sound at the center of the core is 



B4r) 



2mc 



where 



n 1 ; '^"o 

ilo — —UnkuJ — -UJ. 

2 2d 



(26) 



(27) 



It can be shown that near the vortex core, V x By oc 
and V X cx r^. The expression under the integral in 
Eq. ((24)) is therefore proportional to near the center 
of the core and to rKi{r/X) at r > f . Thus, the integral 
in Eq. (|24|) falls off very rapidly inside the core and can 
be neglected. 

The case of fcA > 1 is rather involved as it requires 
explicit knowledge of the structure of the vortex core. 
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For fcA < 1 Eq. HT]) provides that = -(2mc/e)f2 in 
all regions of space, so that A^i 0. In this limit the 
Meissner field Bm and the fields due to images can be 
neglected. The total interaction energy per unit length 
of the vortex acquires the simplest form atK = A/^::^I: 



A^2 



2tt 



fcA 

K 



In K , 



(28) 



where L is the dimension of the slab in the z-direction. 

If one excludes small contribution from the vortex core 
in Eq. ((20|) . then the integration by parts yields 



A2 „ , , 1 

— — = — (bB^x (VxB^,)-ds- — 
h on I 47r 



d^rH B, 



(29) 



This approximation is good if A and d are large compared 
to the coherence length ^. Then the vortex energy per 
unit length is 



AE 

~17 



{Anxy- 



■In K 



47r 



(30) 



The first term in this expression is the self-energy of the 
vortex, while the second term is the energy of the inter- 
action of the flux quantum with the external field. 

The condition for the nuclcation of the vortex, AJF2 + 
AE = Q, yields 



where 



e = 1 - 



H 
W 



cl 



(31) 



(32) 



and Hci = $0 In K/(47rA^) is the first critical field that 
follows from Eq. ^ ai AE = 0. Substituting Eq. ^ 
into Eq. ([51]) . one finds the conditions on the frequency 
/ and amplitude mq of the sound needed to nucleate a 
vortex in the geometry shown in Fig. [T] 



/ = 



2d 



Uo 



d_ 



mv 



(33) 



While the last formula was derived under the conditions 
TT^ < ttA <C d, our numerical analysis shows that it 
holds even for d ^ ttA at k ;3> 1 and is true by order 
of magnitude for k ^ 1. In this case the expression 
for Hci carries the signature of the surface barrier [l^ : 
Hci = /3$o lnK/(47rA^), where 



/3 



1 - 2sinh(d/2A)[sinh(d/A)]-i 



(34) 



For the speed of the transverse sound v ^ ix lO^cm/s, 
in a slab of thickness d ~ ttA ~ 6 x 10~^cm and k ~ 2, 



with H within one percent of Hd , one gets from Eq. p3|) 
/ ~ 3GIIz and uq ~ 0.2nm. These are accessible values 
of frequency and amplitude of ultrasound. 

As r2 changes its sign every half a period of the sound, 
vortices are periodically attracted and repelled by the 
standing acoustic wave in Fig. [TJ Periodic entering and 
expulsion of vortices should result in the elevated atten- 
uation of the ultrasound and in the ac voltage across the 
slab at the sound frequency. In a different experiment 
one can assist vortices to enter or exit the superconduc- 
tor with the help of the surface acoustic waves (SAW). 
Like in the problem with a slab, local rotation of the 
crystal produced by the SAW may assist nucleation of 
the vortex at the field just below Hci- 

In conclusion, we have demonstrated that nucleation 
of a vortex in a superconductor can be assisted by ultra- 
sound. In the presence of a standing sound wave, vor- 
tices will periodically enter and exit the superconductor. 
The required amplitude and frequency of ultrasound are 
within experimental reach. 
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